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1 . INTRODUCTION 

The  problem  considered  in  this  report  concerns  the  relationship 
between  components  of  the  measured  fields  due  to  an  electromagnetic 
pulse  (EMP)  simulator  and  the  possibility  of  extracting  information 
about  the  conductivity  and  permittivity  of  the  ground. 

Radiation  fields  emitted  by  a pulser  can  locally  be  approximated  by 
plane  waves,  and  the  monochromatic  components  of  the  incident  and 
reflected  pulses  are  related  by  Fresnel's  coefficients.  These 
coefficients  involve  the  ground  parameters;  when  they  are  known,  the 
free  field  can  be  computed  from  the  resultant  fields  in  the  air. 
Alternatively,  two  measurements  of  different  fields  for  the  same  pulse 
can  be  compared  for  different  values  of  the  parameters  until  agreement 
is  found. 

Simulators  such  as  the  Transportable  Electromagnetic  Pulse 
Simulator  (TEMPS)  and  the  Army  EMP  Simulator  Operation  (AESOP)  produce 
mainly  horizontally  polarized  fields,  especially  on  the  center  line  (see 
fig.  1) . For  these  fields,  the  reflected  electric  field  is  almost 
180  deg  out  of  phase  with  the  incident  field,  and  the  peak  of  the 
resultant  field  is  determined  mostly  by  the  delay  time, 
returning  quickly  to  zero  after  that.  The  vertical  component  of  the 
resultant  magnetic  field  is  strictly  proportional  to  the  electric  field. 
Simultaneous  measurements  of  these  two  fields  can  be  used  to  determine 
that  the  respective  probes  produce  traces  that  are  in  agreement  with 
each  other.  The  horizontal  component  of  the  magnetic  field  is 
proportional  to  the  difference  (essentially  the  sum  of  the  magnitudes) 
between  the  incident  and  reflected  fields  and  shows  much  more  detail  of 
the  incident  pulse.  This  component  can  be  used  to  determine  the  free 
field  and  then  the  resultant  electric  field,  by  using  the  Fresnel 
coefficient.  When  the  correct  values  of  the  ground  conductivity  and 
dielectric  constant  are  used,  this  procedure  serves  to  verify  both  the 


Figure  1. 


The  TEMPS  coordinate 


system. 


ground  parameters  and  the  agreement  of  the  probes.  Furthermore,  precise 
measurements  of  the  field  would  allow  a determination  of  the  frequency 
dependence  of  the  ground  parameters. 

Two  sets  of  field  measurements  were  used  in  this  work.  One  was 
taken  with  the  sensor  box  on  the  ground;  the  other,  at  a height  of  5 m. 
’’he  quality  of  the  photographs  was  acceptable,  but  much  improvement  is 
needed . 

The  electric  field  was  computed  from  the  radial  magnetic  field,  and 
this  curve  was  compared  to  the  measured  electric  field.  The 
corresponding  spectra  also  were  compared.  The  values  of  the  qround 
parameters  were  varied  to  find  what  effect  they  have  on  the  agreement 
between  both  sets  of  curves,  and  reasonable  values  of  the  parameters 
were  found  that  produced  acceptable  matching  of  the  curves.  Other  input 
parameters,  such  as  the  height  and  the  calibration  constants  of  the 
probe,  were  varied,  but  no  significant  improvement  was  found  to  justify 
such  a change. 

A somewhat  different  way  of  comparing  the  curves  also  is  presented, 
and  it  should  be  tried  in  future  work.  It  could  avoid  some  difficulties 
with  the  numerical  computations. 

Similar  considerations  are  valid  for  "vertical"  polarization.  The 
vertical  component  of  the  resultant  electric  field  and  the  horizontal 
magnetic  field  are  proportional  to  each  other,  and  they  can  be  used  to 
compute  the  radial  component  of  the  electric  field.  Field  mapping  has 
shown  that  the  electric  field  produced  by  the  TEMPS  has  a sizable 
vertical  component  off  the  center  line. 

Once  this  procedure  is  well  established,  it  can  be  used  to  map  the 
ground  parameters  at  a given  site  and  to  observe  their  variation  with 
the  water  content  and  other  local  changes.  More  refined  measurements 
also  could  be  used  to  determine  the  frequency  dependence  of  the  ground 
parameters  for  the  range  of  interest  for  EMP . 

Information  on  numerical  Fourier  transforms  can  be  found,1  and  the 
theory  of  pulses  incident  on  a conductina  qround  has  been 
published.'0  The  TEMPS  field  mapping  is  available,1  and  details  about 
the  digitization  of  the  traces  will  be  the  subject  of  a future  report  by 
Thomas  V . Noon . 

1 Alfred  G.  Brandstein  and  Eqon  Marx,  Numerical  Fourier  Transform, 
Harry  Diamond  Laboratories  TR-1748  (September  1976)  . 

^ Fgon  Marx,  Reflected  and  Transmitted  Fields  for  a Plane-wave  Pulse 
Incident  on  Conducting  Ground,  Harry  Diamond  Laboratories  TR-1740  (April 
1976)  . 

E.  Patrick  and  S.  Soo  Hoo,  Transiortable  Electromagnetic  Pulse 
Simulator  (TEMPS)  Preliminary  Field  Mapping  Report,  Harry  Diamond 
Laboratories  TM-74-15  (October  1974) . 
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FIELD  COMPONENTS 


An  EMP  simulator,  such  as  the  TEMPS  or  AESOP,  produces  a pulse  that 
on  the  center  line  is  horizontally  polarized  to  a good  degree  of 
approximation.  Locally,  the  fields  correspond  to  those  of  a plane  wave, 
and  Fresnel's  formulas  can  be  used  to  determine  the  reflected  wave.  The 
ground  is  assumed  to  be  plane  and  homogeneous. 

To  determine  the  reflected  field,  the  incident  field  is  decompose'.' 
into  its  Fourier  components, 


1 f lWt 

= — / E e dw 

2tt  J a) 


Eft) 

then  the  reflected  field  is  obtained  from 


/°° 

_ _ , , iwt, 

EA(w,e  dw  ' 


sin  iJj  - 

[K  ‘ 1 (o/EoW) 

I - COS2lJjJ 

\h 

sin  t + 

[ ~ i(a/e0u) 

- COS^ll 

h 

The  total  field  in  the  air  at  height  h is  then  given  by 

Ey<t>  - Eh(t>  t E'  (t  - tD)  , 

where  the  time  delay  is 


(1) 


E 1 ft)  = cft  / E..R^(w)e~“''dw  , (2) 

where  the  Fresnel  coefficient  for  horizontal  polarization  is 


(3) 


(4) 


t = 2h  sin  i p/c  . 


(5) 


The  magnetic  field  of  each  pulse  is  perpendicular  to  the  electric  field, 
and  the  vertical  and  horizontal  components  are 


«x(t)  = Zo'  [Eh(t)  " Eh  (fc  ' tD)]sin  * ' 


H (t) 
z 


- z;’  [vt!  * Eh  (‘  - *■>)] 


cos  >Jj  , 


(6) 


(7) 


where 


Z = (u  /e  )**  . (8) 

Thus,  the  vertical  component  of  the  magnetic  field  is  proportional  to 
the  electric  field.  For  a perfect  conductor,  R^  = -1,  and  these 
components  follow  the  free  field  until  time  t and  then  rapidly  return 
to  zero,  due  to  the  cancellation  of  the  fields  of  the  incident  and 
reflected  pulses.  On  the  other  hand,  the  horizontal  component  of  the 
magnetic  field  is  proportional  to  the  delayed  sum  of  the  two  pulses  and 
approximates  the  shape  of  the  free  field  when  the  time  delay  is  small. 

The  Fourier  transform  of  the  horizontal  component  of  the  magnetic 
field  is 

H = Z"1  (l  - p e"1WtD)  sin  ^ E , (9) 

xio  o V h 7 uj 


and  the  free  field  can  be  obtained  from  the  measured  H . Furthermore, 

x 


E = (l  + R,e“1U)tD)E  (10) 

yw  ' n ' w 

also  can  be  computed,  and  its  Fourier  transform  should  reproduce  the 
measured  electric  field. 


3.  FIELD  MEASUREMENTS 

Two  sets  of  measurements  were  taken  for  this  report.  They 
consisted  of  simultaneous  measurements  of  E and  H with  different  faces 
of  the  same  Stanford  Research  Institute  dome- type  Censor  box.  Different 
sweep  speeds  were  recorded  during  separate  shots;  that  is,  only  two 
pictures  were  obtained  from  each  shot.  The  measurements  were  taken  c.n 
the  center  line,  at  100  m from  the  simulator.  The  first  set  was 
obtained  from  a sensor  on  the  ground  with  AESOP;  the  second  set,  at  a 
height  of  5 m with  the  TEMPS  pulser  in  the  AESOP  frame. 

The  time-amplitude  traces  are  shown  in  figures  2 to  5.  The 
composite  digitized  traces  of  the  magnetic  field  and  the  frequency 
spectra  are  shown  in  figures  6 to  9 . 
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Figure  6.  Composite  time-amplitude  trace  of  horizontal  magnetic  field  on 
ground. 


Figure  7.  Frequency  spectrum  of  horizontal  magneto.  f i . id  on  ground. 


4. 


NUMERICAL  PROCEDURES 


For  use  as  input  to  the  computer  program,  the  time-amplitude  traces 
have  to  be  digitized.  Each  trace  is  digitized  separately,  and  a 
computer  proqram,  TRANS,  developed  by  Thomas  V.  Noon,  at  the  Harry 
Diamond  Laboratories  (HDL) , gives  a single  set  of  coordinates  for  the 
entire  curve. 

Several  problems  are  associated  with  this  procedure.  It  is  assumed 
that  pulses  in  a set  are  identical,  but  they  are  not  really;  normally, 
deviations  are  small  enough  so  that  no  severe  errors  are  introduced  in 
this  manner.  A point  has  to  be  identified  in  each  picture  that 
corresponds  to  the  end  of  the  previous  trace;  this  identification  is 
sometimes  difficult,  but  an  error  that  was  inadvertently  made  was 
discovered  through  discrepancies  in  the  two  curves  for  the  electric 
field.  The  ramp  that  precedes  the  actual  pulse  is  difficult  to  obtain 
for  the  magnetic  field  and  essentially  impossible  for  the  electric 
field;  this  ramp  was  digitized  from  the  third  picture  in  the  first  set, 
and  it  was  assumed  to  be  the  same  for  the  second  set.  Another  problem 
related  to  the  ramp  is  the  lack  of  a precise  time  when  the  pulse  starts ; 
when  comparing  the  two  pulses  in  the  output,  one  pulse  had  to  be  shifted 
until  the  leading  edges  coincided.  More  generally,  the  quality  of  the 
pictures  obtained  from  the  oscilloscopes  is  such  that  there  is 
considerable  uncertainty  in  the  definition  of  the  points  on  the 
digitizer . 

Once  an  acceptable  composite  trace  of  the  magnetic  field  was 
obtained,  the  points  were  linearly  interpolated  to  obtain  an  equispaced 
set  of  points  suitable  for  a Fast  Fourier  Transform  (FFT) . 
Then  equations  (9)  and  (10)  were  used  to  obtain  the  Fourier  transform  of 
the  electric  field,  and  an  inverse  FFT  produced  a computed  trace  for  the 
electric  field,  which  was  plotted  with  the  measured  field.  Since  this 
field  goes  quickly  back  to  zero,  only  the  initial  portion  of  the  trace 
was  considered  of  interest,  and  most  of  the  resulting  ramp  also  was 
eliminated  from  the  plots.  Alternatively,  the  absolute  values  of  the 
Fourier  transforms  of  the  computed  and  measured  electric  fields  were 
plotted  on  a semiloq  plot  either  for  the  full  range  given  by  the  FFT  or 
for  only  the  initial  part. 

A decision  that  has  to  be  made  is  what  the  number  of  interpolated 
points  will  be.  The  number  8192  was  chosen,  because  it  samples  the 
curve  adequately  in  the  reqion  with  the  largest  density  of  points  when 
the  final  time  is  chosen  about  the  end  of  the  magnetic  field  pulse,  at 
7.5  or  10  ns.  Doublinq  or  quadrupling  this  number  of  points  did  not 
introduce  any  significant  change  in  the  resulting  trace. 


A serious  problem  that  arose  was  a large  oscillation  superimposed 
on  the  computed  trace  of  the  electric  field.  With  the  help  of 
William  T.  Wyatt  of  HDL,  the  reason  for  this  oscillation  was  found  to  be 
t'.a  division  by  the  factor 

i - 

to  obtain  the  free  field.  The  value  of  R,  is  close  to  -1,  and 


1 + e-iwtD 


vanishes  whenever 


wtD  = kn  , k = 1,  3,  5,  . . . , (11) 

which  causes  a peak  in  the  computed  E for  frequencies  close  to  this 
value.  The  numerator,  Hxt0,  also  should  vanish  for  these  frequencies, 
but  recording,  digitization,  and  computation  errors  account  for  the 
appearance  of  a large  peak  in  E . Similar  computations  done  for 
analytic  curves  (such  as  a double  exponential),  displaced  in  time,  show 
basically  the  same  behavior. 

For  the  box  on  the  ground,  the  time  delay  is  0.85  ns,  and  the  first 
frequency  at  which  this  peak  occurs  is  589  MHz.  With  8192  points  and  a 
time  of  7.5  jjs,  the  significant  frequency  limit  for  an  EFT  is  546  MHz, 
and  this  problem  barely  shows  up.  On  the  other  hand,  for  a 5-m  height, 
the  time  delay  is  4.9  ns,  and  the  frequency  of  the  first  peak  is  only 
101  MHz. 

The  solution  that  was  used  for  this  difficulty  was  a filter 
function  that  multiplied  E . Several  were  tried,  and  the  function 
finally  chosen  was 

F((d)  = | sin  cDtp/wtpl1"1  , (12) 


which  is  1 for  small  <d  and  vanishes  if  u>  is  given  by  equation  (11)  . The 
value  n = 2 was  used  for  the  box  on  the  qround;  n = 1,  for  the  5-m 
height.  If  this  filter  is  used,  the  frequency  spectrumaround  the  first 
zero  and  after  that  is  virtually  ignored,  due  to  the  o dependence. 

Another  method  that  could  be  tried,  especially  for  a height  of  5 m, 
would  be  based  on  the  clear  time  of  4.9  us  to  obtain  the  corresponding 
reflected  field,  subtract  it  from  the  measured  trace,  and  proceed  in 
this  way  to  obtain  the  free  field.  This  method  was  not  tried. 
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It  is  possible  also,  following  a procedure  used  by  Janis  Klebers  of 
HDL,  to  use  equations  (4)  and  (6)  to  separate  the  incident  and  reflected 
fields , 

E ( t)  = E (t)/2  + Z H (t)  / (2  sin  i|> ) , (13) 
h y ox 


E'  (t  - t \ = -E  (t)/2  + Z H ( t)  / ( 2 sin  il)  . (14) 
hV  D / y ox 

Then  the  matching  of  the  Fresnel  coefficient  can  be  done  from  the 
Fourier  transforms  of  E^  and  E'.  One  difficulty  with  this  approach  is 
the  need  to  know  a common  initial  time  for  the  traces,  which  knowledge 
is  difficult  to  get  because  the  ramp  is  not  recorded  for  the  electric 
field. 


5.  VARIATION  OF  PARAMETERS 


A number  of  different  parameters  can  be  varied  in  this  computer 
program.  The  main  concern  of  this  work  related  to  the  conductivity  of 
the  ground,  o,  and  the  dielectric  constant,  k.  The  height  of  the  sensor 
also  was  varied,  as  were  the  probe  calibration  constants,  but  the 
measured  values  were  found  to  be  adequate.  The  angle  of  incidence  was 
determined  from  the  geometry  of  the  experiment. 


A good  coincidence  of  the  computed  and  measured  electric  fields  for 
the  sensor  at  a height  of  0.65  m aboveground  was  obtained 
with  o = 20  mmho/m  and  ►:  = 25  (fig.  10,  11).  Figure  12  shows  the  full 
range  of  the  frequency  spectrum  that  is  obtained  from  an  8192-point  FFT . 
There  is  good  agreement  in  the  shape  of  the  Fourier  transforms 
throughout  the  range,  even  though  the  high-frequency  end  of  the  computed 
transform  is  incorrectly  enhanced  by  the  division  by 


+ V 


-iu)tD 


There  is  an  unexplained  dip  in  the  transform  of  the  measured  electric 
field  at  about  280  MHz,  where  the  computed  spectrum  has  a peak. 
Agreement  in  the  traces  is  improved  by  filtering  out  this  frequency,  by 
using 

F(io)  = (sin  uti/uti)  (sin  wt  , (15) 


where  t ; = 1.732  ■ 10~  s,  very  nearly  twice  the  time  delay. 
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Figure  12.  Full  range  of  frequency  spectrum  as  provided  by  Fast  Fourier 
Transform,  for  o = 20  mmho/m  and  k = 25. 


Figures  13  to  20  show  the  effects  of  varying  1 and  ■ to  values  at 
the  extremes  that  can  be  expected  for  the  ground.  The  conductivity  is 
changed  from  0.5  to  100  mmho/m;  the  dielectric  constant,  from  1 to  100. 
The  effects  of  these  changes  are  quite  noticeable;  the  conductivity 
affects  mainly  the  late  times  and  loi;  frequencies,  while  the  dielectric 
constant  changes  mainly  the  peak  of  the  time-amplitude  trace  and 
somewhat  higher  frequencies. 


Figures  21  and  22  show  the  effects  of  the  change  in  the  assumed 
height  of  the  sensor  aboveground,  which  through  the  time  delay  affects 
mainly  the  peak  in  the  time- amplitude  trace. 

Figure  23  shows  the  computed  trace  that  is  obtained  by  this  method 
when  no  filter  is  used,  and  figure  24  shows  the  trace  for  a 4096-point 
transform  with  no  filter.  The  latter  is  less  accurate,  but  is  not 
affected  by  the  "resonance,"  because  it  does  not  go  as  high  in 
frequency.  Figure  25  shows  the  effect  of  the  filter  in  equation  (15). 
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Figure  13.  Time-amplitude  traces  for  o - 0.5  mmho/m 
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Figure  14.  Frequency  spectra  for 
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Figure  15.  Time-amplitude  traces  for  o = 100  mmho/m 
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Figure  16.  Frequency  spectra  for  o = 100  mmho/m. 
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Figure  18.  Frequency  spectra  for  ■ = 2o  mmno  m,  • 


Figure  19.  Time-amplitude  traces  for  o = 20  mmho/m,  k 
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Figure  20.  Frequency  spectra  for  o = 20  mmho/m,  k 
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Figure  21.  Time-amplitude  traces  for  o ^ 20  mmho/m,  < = 25,  assuming 
that  height  of  sensor  aboveground  is  0.55  m. 
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Figure  22.  Time-amplitude  traces  for  o - 20  mmho/m,  ► 25,  height 

of  sensor  aboveground  0.75  m. 
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Figure 


'3.  Time-amplitude  traces  for  a = 20  mmho/m,  < = 25,  showing 
spikes  when  filter  is  not  used  with  8192-point  transform. 
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Figure  24.  Time-amplitude  traces  for  1 20  mmho/m, 

filter  and  with  4096-point  transform. 
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Figure  25.  Time-amplitude  traces  for  o = 20  mmno/m,  k = 25,  with 

filter  that  corresponds  to  dip  in  frequency  spectrum  at 
280  MHz. 


Figures  26  to 
5 m aboveground , 
resonance  occurs 
agreement  at  low 
ground . 


28  show  the  same  graph  for  the  measurements  taken  at 
with  the  values  of  o = 30  mroho/m  and  - 25.  The 
at  a correspondingly  lower  frequency,  and  the  genei al 
frequencies  is  not  so  good  as  that  for  the  box  on  the 


Figure  26.  Time-amplitude  traces  for  a = 30  mmho/m,  k = 25,  when 

measurements  are  taken  with  sensors  at  5 m aboveground. 
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Figure  28.  Full  range  of  spectra  for  ' = 30  mmho/m,  ► 25,  and 
height  of  sensor  aboveground  = 1 

6.  CONCLUSIONS 

A good  agreement  was  obtained  for  the  time-amplitude  traces  and  the 
frequency  spectra  for  the  chosen  best  values  of  the  parameters;  however, 
considerable  improvement  is  needed. 

For  both  sets  of  measurements,  on  the  ground  and  at  5 m,  the  peaks 
in  the  electric  field  could  be( matched  quite  well.  The  agreement  was 
inferior  in  the  following  oscillations,  to  improve  again  at  later  times 
before  the  field  essentially  went  back  to  zero.  There  was  no 
significant  qualitative  difference  in  the  matching  for  the  two  sets. 

The  frequency  spectra  showed  good  agreement  also  at  lower 
frequencies,  up  to  200  MHz  for  the  box  on  the  ground  and  up  to  50  MHz 
for  the  5-m  height.  The  wide  divergence  at  higher  frequencies  was  found 
to  be  due  to  the  limited  accuracy  of  numerical  computations. 
Surprisingly,  the  general  shapes  of  the  curves  aqree  quite  well  uf  to 
the  limit  of  550  MHz  for  the  measurements  on  the  ground,  except  that  s 
dip  in  the  spectrum  of  the  electric  field  at  about  280  MHz  was  not 
explained.  The  quality  of  the  agreement  is  significantly  better  for  the 
measurements  taken  on  the  ground  than  for  those  taken  at  5 m. 


A number  of  factors  that  limited  the  accuracy  of  these  comparisons 
can  be  improved  without  much  difficulty.  The  different  sweep  speeds  for 
the  oscilloscopes  can  be  taken  simultaneously  when  the  instrumentation 
van  is  available.  The  quality  of  the  photographs  can  be  improved  to 
show  more  detail  in  rapidly  varying  portions  of  the  curve  and  less 
intensity  for  the  slowly  varying  parts.  The  amplifications  can  be 
increased  to  show  more  detail  at  late  times.  Time  delays  can  be  used  in 
the  triggering  circuit  to  obtain  a picture  with  a high  sweep  speed  at 
later  times,  to  obtain  more  detail  where  necessary.  Other  numerical 
techniques  can  be  tried  to  avoid  some  of  the  difficulties  with  false 
resonances. 

Comparison  of  the  measured  and  "computed"  fields  shows  that  the 
probes  for  the  electric  and  magnetic  fields  agree  quite  well  and  are 
probably  accurate  in  an  absolute  sense  within  the  range  of  interest. 
Also,  the  calibration  constants  for  the  probes  were  used  without 
modifications,  so  at  least  their  ratio  is  correct.  The  assumptions  that 
the  pulse  was  a plane  wave  over  a homogeneous  plane  infinite  ground  in 
the  radiation  zone  for  the  field  also  did  not  introduce  any  large 
errors . 

That  the  wave  comes  from  different  parts  of  the  simulator  and  not 
from  a point  source  at  a large  distance  might  explain  some  of  the 
discrepancies  observed  after  the  first  few  peaks.  Also,  the  effect  of 
the  metal  box  on  the  measured  field  was  not  taken  into  account. 

In  summary,  this  method  provides  useful  information  about  the 
sensors  and  the  ground  constants,  and  it  can  be  extended  to  fields  off 
the  center  line  and  other  combinations  of  components. 
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